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Superconductor/normal conductor/superconductor (SNS) Josephson junctions with highly trans- 
parent interfaces are predicted to show significant deviations from sinusoidal supercurrent-phase 
relationships (CPR) at low temperatures. We investigate experimentally the CPR of a ballistic 
Nb/InAs(2DES)/Nb junction in the temperature range from 1.3 K to 9 K using a modified Rifkin- 
Deaver method. The CPR is obtained from the inductance of the phase-biased junction. Transport 
measurements complement the investigation. At low temperatures, substantial deviations of the 
CPR from conventional tunnel-junction behavior have been observed. A theoretical model yielding 
good agreement to the data is presented. 

PACS numbers: 74.45.-l-c, 74.50.-|-r 



I. INTRODUCTION 

A fundamental item in the theoretical modeling of a 
Josephson junction is the dependency of the supercurrent 
Is flowing through it on the phase difference ip across the 
junction. Golubov, Kupriyanov and Il'ichev published a 
recent review on the subject.^ As this relation is neces- 
sarily 27r-periodic and odd,^ it can be expressed as the 
following Fourier series (with the critical current Ic)- 



(1) 



/s(<^) is known as (super-)current-phase relationship or 
CPR, the dimensionless term f{(p) as normalized CPR. 
For vanishing transparency (e.g. a tunnel junction), only 
Ii is relevant, reducing Eq. ^ to the well-known dc 
Josephsoni equation 



Isi'p) = Ii sin ip = Ic sin ip. 



(2) 



Junctions with direct (i.e. non-tunnel) conductivity, such 
as the Nb/InAs(2DES)/Nb weak Hnks used in this work, 
are predicted to show more complex behavior and a non- 
sinusoidal CPR due to higher order processes of charge 
transport such as multiple Andreev reflect ions Indi- 
cations of a significant non-harmonic term in our junc- 
tions have been found in microwave measurements.^ We 
are interested in experimental access to the CPR as a way 
to test and improve the theoretical models describing our 
junctions, since the conduction mechanisms in these de- 
vices are still not fully understood. The measurements 



deliver not only the dependency f{(p), but also Ic, thus 
offering an independent comparison to transport mea- 
surements. 

II. THEORY 



The first extensive theoretical works concerning high 
transparency superconducting weak links were published 
by Kulik and Omel'yanchuk for a short quasiclassical 
point contact in the dirty limit^ I ^ i^o (with mean free 
path I and coherence length i^o) and for the clean limit^ 
/ 3> ^Oj predicting a CPR f{ip) changing from a sinusoidal 
curve to a saw tooth shape for high transparencies and 
low temperatures. These idealized model systems serve 
to understand the basic mechanisms leading to a devia- 
tion of the CPR from the well-known Josephson relation 
(Eq. 0), which we encounter also in complex real-world 
junctions, though mixed with secondary processes. Me- 
chanically controlled break-junctions* come closest to the 
assumptions of the theory of Kulik and Omel'yanchuk. 
CPR measurements by a flux-detecting method in this 
system^ agreed with the predicted changes in the posi- 
tion of the maximum and curve shape. In contrast to 
these systems, our SNS junctions exhibit a finite length, 
some scattering of electrons in the 2DES, and an area of 
induced superconductivity around the electrodes caused 
by the proximity effect. To take these properties into 
account and to describe the non-sinusoidal CPR in our 
experiments, we used the scattering matrix formalismifi 
to derive a formula for the Josephson current through a 
double-barrier structureiii 
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where Rq = h/e^ is the resistance quantum, A is the coef- 
ficient of the Andreev reflection at the InAs(2DEG)/InAs 
interface and Sij are the elements of the scattering ma- 
trix for normal reflection in a symmetric double bar- 
rier junctionl^ Sn = |r| + \r\\t\'^pl/{l + \r\'^pl), S12 = 
+ kPp^i)- Here |rp, |tp are the reflection and 
transmission probabilities of the left and right-hand bar- 
rier, pn — eiq){ikFnCL — tOma/hvpn), where uj„i — {'^rn + 
Vj'KksT are Matsubara frequencies, kpn and vpn are 
components normal to the barriers of the wave vector 
and Fermi velocity, respectively, of the n-th transverse 
mode to the barriers and a is the distance between them. 
Following the BTK approachji^ the reflection and trans- 
mission coefficients can be written in the form 



(l-(n/iV)2)(r; + 1)2/477 + ^2 w 

where Z is the dimensionless potential barrier strength 
and 77 = vps/vpn is the Fermi velocity mismatch. 

In order to calculate the Josephson current from Eq. O 
one must determine the Andreev reflection coefficient A 
for the InAs(2DEG)/InAs interface. Since superconduc- 
tivity is induced in InAs by the proximity effect, the co- 
efficient A may be written aa— A = iF/{l + G), where 
F and G are Green functions in the inversion layer of 
InAs. Due to the low electron density of the inversion 
layer the suppression of the pair potential in Nb can be 
neglected and F, G can be expressed by the McMillan 
equationsi^ G = lo/Vuj'^ + F = ^/y/uj'^ + $2, and 
$ = A/(l -I- 7s\/'i2 -I- A2), where 7s is a dimensionless 
parameter characterizing the transparency between Nb 
and InAs, A = A/irkBTc, lu = uj/TrkBTc, A is the su- 
perconducting energy gap of bulk Nb, and Tc its critical 
temperature. Both the critical temperature and the en- 
ergy gap of the Nb can be suppressed near the interface 
because of disorder— but A « 0.6 remains constant— 
Since Tq can be determined from temperature measure- 
ment, there are only two free fitting parameters in the 
model, the carrier density rig and 75 (we assume Z = 0, 
i.e. no real barrier). The Fermi velocity vp and Fermi 
wave vector kp are calculated for a given value of n^. 
Since the normal resistance depends on as well, the 
Us obtained from the fit can be verified comparing the 
theoretical and experimental value of the resistance. 
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FIG. 1: Temperature dependence of the critical current. Cir- 
cles are experimental data, the solid line presents a fit using 
the model described in the text. The inset shows a schematic 
cross section of the Josephson junction. 

lapping sample geometry (see inset of Fig.^, where we 
can set the electrode separation a with nm accuracy using 
anodic oxidation to grow the insulating Nb205 interlayer. 
At a typical electron density of ris = 1.2 x 10^^ cm~2 
and mobility of /i « 10"*^ cm2/Vs, the mean free path of 
I = 240 nm is much longer than the electrode separation 
a = 40 nm. From the dependence of the critical current 
on the channel length, we estimate the coherence length 
£,N ~ 145 nm at 1.8 K. Thus we have a ballistic (a <C I) 
superconducting weak link in the short limit (a <^ £,n)- 
More details have been published by Chrestin et al— 

The use of the native 2DES on p-type InAs is not with- 
out difficulties, as it forms on every surface of the crys- 
tal. Thus it offers parallel conduction paths, resulting 
in significant bypass currents around the junction, which 
reduce the normal resistance Rn. 



IV. CPR MEASUREMENT TECHNIQUE 



III. SAMPLE PREPARATION 

The superconductor used in our SNS junctions is a Nb 
thin film, the normal conductor is the two-dimensional 
electron system (2DES) that forms as naturally occur- 
ring inversion layer at the surface of bulk p-type InAs. 
The InAs surface is cleaned in situ using low-energy Ar 
etching prior to the deposition of a 100 nm thick Nb film 
by magnetron sputtering, yielding the highly transparent 
interfaces required to observe deviations from sinusoidal 
CPR. The structures are patterned using standard opti- 
cal and electron-beam lithography. We are using an over- 



The CPR of the Josephson junction is determined by 
an inductive rf frequency method requiring no galvanic 
contacts to the sample, thus reducing noise. This tech- 
nique is based on the work of Rifkin and Deaver JiS, The 
junction to be investigated is incorporated into a super- 
conducting loop. This loop forms an rf SQUID of induc- 
tance L, which is coupled inductively to a high quality 
tank circuit in resonance. The phase difference if across 
the junction can be biased by an external magnetic flux 
<I>e via a dc current Idc- Changes in the impedance of 
the coupled system are measured and can be used to re- 
construct the CPR. Further details are given in Ref. Q. 
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FIG. 2: Measured phase shift at various temperatures. 



This method requires low critical currents Ic of the junc- 
tion and low inductances L of the SQUID washer, as the 
SQUID enters the hysteretic regime when Pf'{(p) > 1, 
where p = 2ti Lie I 'J'o is the normalized critical current 
and f'i^p) = df{ip)/dip. In the hysteretic mode the in- 
ternal magnetic flux as a function of the external one 
becomes multivalued, preventing us to reconstruct the 
CPR in the complete phase range (0, 2tt). The challenge 
is to reduce Ic and L while maintaining excellent inter- 
face transparency, large IqRn products, and sufficient 
coupling to the tank circuit. We achieve this by reducing 
the geometric dimensions of the junction and the loop 
and by complex washer designs. The presented junction 
is connected to six SQUID washers in parallel with an in- 
tegrated flux transformer, for an inductance of L ^ 17 pH 
and optimized coupling to the tank circuit. A similar 
transformer is described and depicted in Ref. 1. 



V. RESULTS 

The experiments are performed on an overlapping 
Josephson junction of width w = 1 /im and electrode 
separation a = 40 nm. Transport measurements at 1.8 K 
result in Ic = 3.8 /lA, Rn = 58 fl, IcRn = 220 /xV, 
and show clearly developed subharmonic gap structures 
(SGS) in the differential resistance of the junction. This 
indicates high interface quality and transparency of the 
junctions. 

Figure 121 shows the recorded phase shift a in the tank 
circuit as a function of the applied quasi-dc current /dc 
ramping the external flux $e through the SQUID loop. 
This signal was averaged for 20 periods to reduce noise 
and then used to reconstruct the CPR depicted in Fig.O 
For decreasing temperatures, gradual deviation of the 
CPR from conventional sinusoidal behavior towards a 
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FIG. 3: CPR (red) at temperatures between 1.3 K and 6.0 K, 
reconstructed from the data in Fig.|5| The solid lines are a fit 
according to the model described in the text. For comparison, 
a pure sine curve is included. 



sawtooth-shaped curve is observed. The position of max- 
imum current is shifted to O.GStt at 1.3 K. At the same 
time, the critical current increases substantially, so that 
at some temperature the SQUID enters the hysteretic 
regime as l3f'{(p) approaches unity. In that case we are 
no longer able to reconstruct the CPR at lower temper- 
atures. 

The Ic — T dependence was fitted by the model de- 
scribed above using the least squares method with two 
fitting parameters Ug and jb ■ The other parameters used 
in the calculation were determined independently as crit- 
ical temperature Tc = 6.5 K, width w = 1.0 /zm, elec- 
trode separation a = 40 nm, interface parameter Z — 0, 
and Fermi velocity mismatch rj — 0.93. The fitting pro- 
cedure yields the parameters Ug — 3.75 x 10^^ cm^^ and 
7b — 2.4 and good agreement between theoretical and 
experimental data (see Figs. ^ and ^ . Note that the 
model fits the current-phase characteristic quite well for 
the parameters determined from the temperature depen- 
dence of the critical current. Only the curve taken at 
T « 2.0 K exhibits some deviation from the theoreti- 
cal one. However, the temperature dependence of Ic{T) 
is rather steep at low temperatures and a small differ- 
ence between the temperature of the thermometer and 
the sample could be responsible. From the above param- 
eters we calculate the junction resistance of Rn = 83 Vl. 
Transport measurements yield a slightly lower value of 
58 f2, but as we measure a parallel connection between 
substrate and junction, a higher junction resistance is to 
be expected. A realistic substrate resistance of 180 VI 
produces the measured value of 58 O. The interface pa- 
rameters indicate a highly transparent junction. The 
assumed carrier density is about three times the result 
TT-s = 1.2 X 10^^ cm~^ from Shubnikov-de Haas measure- 
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ments on comparable samples. Nevertheless, Ug as de- 
termined from our model is consistent with the smaller 
experimental value of the junction resistance. A possi- 
ble explanation for both effects is an effective widening of 
the contact due to the proximity effect in the surrounding 
2DES. 



junctions and the measurements on Josephson field ef- 
fect transistors presented in Ref. 0. A model yielding 
good quantitative agreement to the results is presented. 
Transport measurements support the results gained by 
phase sensitive measurements. 



VI. CONCLUSION 

We have successfully measured the current-phase re- 
lationship (CPR) of a Nb/InAs(2DES)/Nb Josephson 
junction in dependence of temperature. At low temper- 
atures, substantial deviations of the CPR from a sinu- 
soidal behavior towards a saw tooth shape are observed. 
This is in qualitative agreement with predictions of the 
Kulik-Omel'yanchuk theory for highly transparent SNS 
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